From thermally activated to viscosity controlled fracture of biopolymer hydrogels 
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We report on rate-dependent fracture energy measurements over three decades of steady crack ve- 
locities in alginate and gelatin hydrogels. We evidence that, irrespective of gel thermo-reversibility, 
thermally activated "unzipping" of the non-covalent cross-link zones results in slow crack propaga- 
tion, prevaling against the toughening effect of viscous solvent drag during chain pull-out, which 
becomes efficient above a few mm.s"'^. We extend a previous model [Baumberger et al. Nature 
Materials, 5, 552 (2006)] to account for both mechanisms, and estimate the microscopic unzipping 
rates. 

PACS numbers: 62.20.mm, 83.80.Kn 



A variety of biopolymers (e.g. proteins and polysac- 
charides) can self-assemble into physically cross-linked 
networks in aqueous solutions. The resulting hydrogels 
are usually bio-compatible and show mechanical prop- 
erties that mimic those of extracellular matrices. They 
are therefore good candidates as scaffolds for in vivo tis- 
sue regeneration, e.g. cartilage (ij. Such applications 
require sufficient mechanical strength to withstand ma- 
nipulations associated with implantation and in vivo exis- 
tence [21] . Surprizingly, understanding the physical mech- 
anisms at work during the ultimate behavior (yield and 
failure) of non-covalently cross-linked hydrogels is still in 
its infancy. In this context, Baumberger et al. [3]] have 
identified a basic dissipative process in the slow frac- 
ture of gelatin which they accounted for in a minimal 
model: gelatin gels belong to a wider class of biopoly- 
mer networks with extended cross-link zones, e.g. H- 
bond stabilized multiple helices, distributed along the 
chains and acting as mechanical fuses which may un- 
zip or unreel under tension forcing overall polymer 
pull-out without chain scission, at the cost of the vis- 
cous dissipation due to the solvent drag on the whole 
chain. This highly efficient mechanism makes these nom- 
inally weak gels amazingly resistant to crack propaga- 
tion even at moderate velocities of 1 mm.s~^ and above. 
In contrast, under quasi-static loading physical gels gen- 
erally creep under stress 0] and ultimately break after 
a strongly stress-dependent, random delay 0]. Though 
this low rate behaviour, especially in thermoreversible 
networks, undoubtly pinpoints the role of thermal acti- 
vation, it remains unclear whether the rate-limiting pro- 
cess in fracture of such intrinsically disordered materials 
is the nucleation of unstable cracks Q or their subsequent 
growthQ. Generally speaking, unraveling nucleation and 
growth effects is a subtle task Q and, when feasible ex- 
perimentally, studying low velocities, steady state crack 
propagation provides more straightforward insight into 
the fracture mechanisms. Dealing with gelatin gels, how- 
ever, reliable measurement of both small crack-tip veloc- 
ities and the associated energy release rates is hindered 
by two intrinsic limitations: (i) Thermal activation pro- 



motes crosslink rearrangements in the bulk as well, hence 
stress relaxation which competes with crack propagation 
to release elastic energy; (ii) Strong pinning of the crack 
front by network inhomogeneities may result in front in- 
stabilities which makes tip position and fracture area 
measurements ambiguous [10|. As reported in this letter 
we circumvent these drawbacks by using alginate gels in 
which cross-linking is ensured by ionic bridges along ex- 
tended zones between polyelectrolyte chains, with bind- 
ing energies intermediate between H-bond and covalent 
ones. Though non-thermoreversible, these gels can be 
termed physical since interchain bounds remain weaker 
than backbone ones. Taking advantage of their negligi- 
ble creep under small stresses ll| and of the absence of 
crack-front instability down to velocities in the 10 ^m.s""'^ 
range, we clearly evidence that at low enough rates, crack 
propagation in these gels is thermally activated; more- 
over, after checking that they also exhibit a viscosity con- 
trolled regime, we extend a previous model Q to account 
for thermally activated, stress-aided cross-link dezipping. 
The outcome is an estimate of the debonding rate which 
lends further support to the prevailing picture of "egg- 
box" binding zones in alginate (l^ . 

Alginates are polysaccharides composed of sequences 
of two sugar residues, referred to as G and M [l3|. Ac- 
cording to the "egg-box" model [l2| , gelation occurs via 
interchain chelation of divalent ions (eggs), here Ca^+, at 
specific sites (boxes) formed by the assembly of two sub- 
sequences made of several contiguous G- units (Fig. [1]). 
The length of a binding GG pair is a = 0.9 nm. 

Sodium alginate was purchased from Kalys (France). 
From supplier's specifications (average molecular weight 
216 kDa and 55% "M" residues) the average contour 
length is A = 550 nm. We use different alginate con- 
centrations c between 0.5 and 2.5 g for 100 ml of solvent 
(deionizcd water, otherwise specified) and a constant ra- 
tio of 10^'^ mol of Ca^+ for 1 g of alginate known to 
minimize gel shrin king (syneresis). Homogeous gel sam- 
ples are obtained [1J| by in situ progressive release of 
Ca^"*" ions from insoluble CaCOs particles subsequent to 
the addition of a slowly hydrolyzing acid ("GDL" from 
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Sigma) quickly mixed with the pregel solution before 
casting it in a mold. The samples — 30 cm long slabs, 
1 cm thick and 3 cm wide, attached to parallel grips — 
are stretched along their widths. Details of the set-up 
have been given previously 0]. All experiments are per- 
formed at T = 295 K. Cracks are initiated by notching 
in the mid-plane. When micron-sized CaCOa particles 
(from Sigma) are used, gelling is slow and homogeneous. 
However, at low c, particle sedimentation occurs lead- 
ing to a toughness gradient as revealed by a tilted crack 
front. With nanoparticles of average size 90 nm (Ameri- 
can Elements, CA), no sedimentation occurs but gelation 
is much faster and, at high c, yields inhomogeneous sam- 
ples as revealed by a non-planar, tortuous crack front. 
We have therefore used a mixture of nano- and micro- 
particles, in such proportions as to yield straight crack 
fronts perpendicular to the slab faces and propagating 
along its mid-plane, a very stringent requirement. 




FIG. 1: Schematic hierachy of length scales in the fracture 
of an alginate gel. (a) Experimental set-up showing the mm- 
scale crack opening (b) Crack-tip region, (c) Egg-box rep- 
resentation of a ionic cross-link, (d) Molecular picture of a 
binding unit. 

Steady crack velocities V are measured by video track- 
ing of the crack tip. In our fixed grip configuration, the 
energy release rate G is imposed. We compute it, ne- 
glecting edge effects, from the total stored elastic energy 
as determined from the force vs. stretching ratio load- 
ing curve measured on an unnotched sample. The small 
strain shear modulus fi = 1.5 kPa for c = 1.5% varies 
approximately as in the studied range. 

Figure displays G{V) curves for alginate gels with 
different c. They share with gelatin ones the principal 
features we associate with a visco-plastic fracture mech- 
anism via scissionless chain pull-out a quasi- linear 
strong growth at high V, extrapolating at F = to a non- 
zero threshold value of order a few J.m~^, with a slope 
increasing linearly with the solvent viscosity rj. This lat- 
ter property is confirmed in a control experiment: using 
a glycerol/ water mixture we induce a 2.5-fold increase of 
?7, resulting is a clear toughening with respect to a pure 
water-based gel. When plotted vs. rjV, however, both Q 
curves collapse over the whole F-range (Fig[2|). 

Beyond these remarkable similarities (irrespective of 
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FIG. 2: Fracture energy vs. tip velocity for different algi- 
nate/water wt. fraction (open symbols). Filled circles are for 
a. c — 1.5 % gel in water/glycerol. Filled squares: same data 
set as a function of velocity scaled by solvent viscosity ratio 
(see text) . Lines are best fits by eq. El 



the thermoreversibility of the gels) , Q {V) curves for algi- 
nate specifically exhibit a systematic round-off at lower 
velocities, which implies no clear threshold for crack 
propagation. This is reasonably ascribable to thermally 
activated unzipping. Since with these systems T cannot 
be varied widely enough to build a significant Arrhenius 
plot, we must rely on modeling to test this hypothesis. 
This leads us to extend our previous model [3;] to ac- 
count for thermal activation. In this description, the 
complex network features are lumped into a few param- 
eters, namely the average chain contour length A, the 
areal density So of chains crossing the fracture plane, 
the size a of a binding unit and the activation barrier 
U. The fracture energy is then computed according to 
the Dugdale-Barenblatt theory |15| assuming a uniform 
stress aup over a small-scale cohesive zone at the crack 
tip. The fracture criterion corresponds to the overall pull- 
out of chains, stretched taut, i.e. to a maximum opening 
A of the tip in the cohesive zone (Fig. [TJj). The fracture 
energy is thus simply: 



g = atipA (1) 

Note that G > I J.m'^ entails aup > 2 MPa > fi. 
This justifies that the chains are almost fully stretched 
in the crack tip vicinity. Accordingly, if v is the frequency 
at which units are released, neglecting re-bonding rate, 
the pull-out velocity reads: = av. 

There is little hope to compute the very shape of 
the cohesive pre-crack due to the strongly non-linear, 
anisotropic elastic field which prevails ahead of the tip, 
where stresses reach values of order a^p ^ ^ [l^. In- 
stead we assume a wedge-shaped tip, which provides a 
simple kinetic relationship : 
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where the fitting parameter a should consistently be ^ 1 
to account for crack blunting [lij . 

Closure of the problem requires relating aup to hence 
to v. The chain tension decreases along the polymer away 
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from the crack edge due to a viscous drag of order rj-d per 
unit of contour length. Let us make a crude estimate of 
the tension /y under which reels yield as: 



So 



(3) 



The stress p stems from the capillary pressure jump 
which tends to suck-in the chain when it is pulled-out 
dry. It is zero when the crack tip opening is wetted by a 
drop of solvent [3] but is p = e^fEo/a for a dry tip, with 
en the free energy of solvation per residue of length a 
(the work done by the capillary force against drawing a 
residue into the gap). 

Following Kramers [l3|, the bond breaking rate in the 
biased binding potential is given by: 



iyoexp[-{U - fYa)/kBT] 



(4) 



where vq is the frequency of attempt to escape over the 
binding barrier. 
Eqs. ©-(m yield : 
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As seen on Fig. [2l the functional form of G{V) in 
eq-® provides an excellent overall fit of experimental 
data for gels of composition c ranging from 0.5 % to 2.5 
%. In order to put this analysis on a more quantitative 
footing, we note that the expression of eq. (O features 
three independent fitting parameters only. Disregarding 
multiplicative constants of order unity in eq. ([6]), and 
setting A and a to their nominal average values, there 
remain four parameters actually unknown : J7, a, V* 
and Go- AQh is determined through an independent ex- 
periment where the crack tip is wetted by a drop of the 
solvent ; for gelatin gels this results in a shifted G{V) 
curve. For alginate gels, we have observed no significant 
effect of wetting the tip with water (with NaCl added to 
equilibrate the Na+ concentration with that of the bulk 
sodium alginate one). We conclude that the polyelec- 
trolyte alginate chains are pull-out in a hydrated state, 
hence AQh — 0. 

The relatively simple eggbox structure of interchain 
binding in alginate prompts us to choose the activation 
barrier height U as an input for data fitting with eq.lO. 
To do so, we compute the electrostatic energy of a Ca^+ 
ion assuming it is involved in purely ionic bonds with its 
nearest neighbours O atoms in the GGoGG cage (Fig. 
[TJi) With an average Ca-0 distance r ~ 2.3A this 
entails U/kBT = ils/r = 12 at 295 K in water where 
the Bjerrum length Ib is 7 A. 
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FIG. 3: Fit (solid curve) of g(y) data using eq.© for (a) 
c = 1.5 % alginate and (b) c = 10 % gelatin. Insets : Residual 
energy 5Q after subtraction of the linear component. Dashed 
line on (a) is the T = K characteristics. 



We now restrict our analysis to a c = 1.5 % gel for 
which we have got the most extensive data set although 
the same qualitative conclusions can be drawn from the 
other concentrations. The corresponding fit is shown on 
Fig. The Griffith, or 'T = K" energy threshold 
is Go — 4.6 J.m~^. A mere extrapolation to 1^ = of 
the apparent linear regime over the experimental veloc- 
ity window underestimates this value. Strictly speaking, 
the purely viscous regime is only reached asymptotically 
for V V* . Here, V* ~ 0.6 m.s^^, out of experimental 
range. One can however check on eq.® that the cor- 
rection to the slope of G{V) due to the activated term 
is negligible for V larger than 1 mm.s~^, i.e. that the 
asymptotic slope can be safely estimated from the low 
velocity (V < V*) data set. The value a ~ 8 obtained 
from this slope is indicative of a strongly blunted tip 

Although the fitting value for V* is strongly sensi- 
tive to approximations in the trial value for U, we claim 
that the order of magnitude of the attempt frequency 

~ 5.5 X 10^ s~^ yields deep insight into the unzipping 
dynamics. As discussed by Evans [3], in liquids, the 
thermal impulses that drive unbinding events are damped 
by viscous coupling to the environment. Accordingly, the 
opening of a GGoGG molecular cage is expected to in- 
duce hydrodynamic flow over a size of order a hence a 
damping rate scaling as kBT/rja? = 5.4 x 10^ s~^ in wa- 
ter at 295 K. This is precisely the order of magnitude 
of i/Q, lending strong support to a simple unzipping sce- 
nario where calcium ions would be released one by one. 



4 



as might be schematically expected from the egg box pic- 
ture of the cross-links. This is, to our best knowledge, 
the first time a dynamical argument is given in favor of 
the egg-box architecture. 

On approaching V*, the activation barrier smoothes 
out and the escape rate is no longer given by eq.(|l]) 
since advection-driven, deterministic debonding events 
becomes increasingly prevalent. In the opposite range, 
as — > 0, re-binding events must become relevant and 
lead to a regime ruled by the slow creep of the cross-links 
themselves, with fy going linearly to zero [5]. Conse- 
quently, eq.dH) becomes unphysical for /y ^ i.e. for 
V £ Vmin — V* exp{—U/kBT). For alginate gels, we 
estimate Vmin — 4/xm.s^^, well below the lower bound 
of our velocity set. This legitimates a posteriori using 
eq. (O over the whole experimental window, which fulfills 
the requirement Vmin <^V <$^V*. 

The previous discussion is based on generic features of 
zipper-like cross-linked gels. It can be therefore expected 
that thermally activated, stress-aided dezipping will be 
all the more efficient in thermoreversible gels. This 
prompts us to reassess, at least qualitatively, the case 
of gelatin. As already mentioned, this requires taking 
special care. First of all, using a c = 10 % gelatin/water 
gel, the crack front instability j9| is pushed down be- 
low Vc = 20/xm.s~^. Moreover, we set-up a procedure 
to correct the fracture energy for stress relaxation. For 
this purpose, we record the crack velocity vs. time in 
response to various crack openings. A twin sample, kept 
unnotched, is submitted to the very same stretching se- 
quence while recording the loading force. The stored elas- 
tic energy is computed assuming that stress relaxation in 
such weak, transient networks occurs via debonding of 
cross-links which eventually rebind at a more favourable 
place, therefore resulting in a drift of the reference state, 
hence in a mere shift of the non-linear stress-strain curve 
[l9| . As shown on figlH this tedious procedure ultimately 
reveals a clear logarithmic behaviour over at least two 
decades, essentially below 1 mm.s~^, which was there- 
fore overlooked in previous studies 

According to the present model, breaking one-by-one 
H-bonds {U ~ 0.1 eV) between peptidic residues dis- 
tant of a = 0.3 nm would result in a logarithmic shift of 
the fracture energy of water-based gels strictly between 
V* m.s-i and Kn„ = V* expi-U/ksT) = 8 cm.s"! 
where we use the conservative value a = 10 for the blunt- 
ing parameter [31 . Both the span of this velocity bracket 
and its absolute location are clearly incompatible with 
the experimental observation. We are therefore led to 
propose that the basic thermal event involves the coop- 
erative debonding of n several subsequent H-bonds along 
a cross-link zone. This we attribute to the strong topo- 
logical constraint imposed by the triple-helix structure 
of the cross- links 20], the unzipping of which requires 
also large-scale unwinding. This cooperative mechanism 
cannot be discriminated from a one-by-one unzipping as 



long as one probes (through Qq) the yield tension fy 
which reads nU/na = U/a. However, the characteristic 
unbinding rates are dramatically affected by the effec- 
tive barrier energy (~ nil) and by the bulkiness of the 
activated unit (~ na). On this respect, a cooperativity 
level of n « 3 would make gelatin gels looking similar to 
alginate ones, everything equal otherwise. 

This analysis confirms the importance of thermally 
activated rate processes in soft matter fracture physics 
[1, [2l[. Our claim that the "subcritical" fracture be- 
haviour is sensitive to distinctive topological features of 
zipper-like cross-linked networks opens the way to a more 
extensive experimental study, taking advantage of the 
wealth of network architecture offered by biopolymer hy- 
drogels. 

We are indebted to C. Caroli for in-depth comments 
and careful reading of the manuscript. 
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